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3 ABSTRACT 

^2 ! We report on a systematic center-to-limb variation in measured helioseismic travel 

times, which must be taken into account for an accurate determination of solar interior 

OO I meridional flows. The systematic variation, found in time-distance helioseismology 

analysis using SDO/HMl and SDO/AIA observations, is different in both travel-time 

magnitude and variation trend for different observables. It is not clear what causes 

. _ this systematic effect. Subtracting the longitude-dependent east- west travel times, ob- 

P^l tained along the equatorial area, from the latitude-dependent north-south travel times, 

Q I obtained along the central meridian area, gives remarkably similar results for different 

-)— > I observables. We suggest this as an effective procedure for removing the systematic 

ci ' center-to-limb variation. The subsurface meridional flows obtained from inversion 

of the corrected travel times are approximately 10 m s^^ slower than those obtained 

^ I without removing the systematic effect. The detected center-to-limb variation may 

'nI" I have important implications in the derivation of meridional flows in the deep interior, 

Q\ I and needs a better understanding. 
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1. Introduction 

An accurate determination of the meridional flow speed in both the solar photosphere and the 
solar interior is crucial to the unde rstanding of solar dynamo and predicting solar cycle variations. 



For example, iDikpati et al.l (|2006|) suggested that the slow-down of the meridional circulation dur- 



ing the solar-cycle maximu m could change the duration of the following minimum and delay the 



onset of the next solar cycle. iHathaway & Rightmirel (|2010|) found that the meridional flow speed 



was substantially faster during the solar minimum of Cycle 23 than during the previous minimum, 
and suggested that this might explain the prolonged minimum of Cycle 23. Using flux-transport 
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dynamo model simulations. iDikpati et al.1 (|2010|) suggested that the prolonged minimum of Cycle 
23 might be due to that the poleward meridional flow extended all the way to the pole in Cycle 
23, unlike in Cycle 22 the flow switched to equator- ward near the latitude of 60°. All these works 
demonstrate that an accurate measurement of the meridional flow speed is very important. The pho- 
to spheric meridional flow speed can be determi ned by tracking certain photospheric features, such 
as magnetic structure s and supergranules (e.g.. iKomm et al.l 1 19931 : iHathaway & Rightmird 12010 : 
Hathaway et al.l 12010 : lOizonl |2004l : ISvanda et al.ll2006h . although it is not quite clear how well 
the motions of these surface features represent the photospheric plasma flows. T he photospheric 



meridional flow sp eed can also b e inferred from Doppler-shif t measurements (e.g.. IHathaway et al 



1 996l : luirichlbo 1 0|) . and recently lUlrichl (|20 1 0|) made extensive comparisons of the meridional flow 
speeds obtained by different methods. The meridional flows in the solar interior are primarily de- 
termined by helioseismology, i.e., b y measuring frequency shifts between poleward and equator - 
ward traveling acoustic waves (e.g., iBraun & Fanlll998l: iKiieger et al.l 120071: iRoth & Stixll2008|) . 
and by use of local he lioseismology techniques, narnely, r i ng-diagram analysis and time-distance 



helio s eismology (e.g..lGiles et al. 



20041 : iGonzalez Hernandez et al.l 



20081). 



19971: IChou & Daill200ll : iHaber et al.ll2002l: IZhao & Kosovichev 



The subsurface meridional flow speeds obtained by the two different local helioseismology 
techniques are in rea sonable agreement for at least the upper 20 Mm of the convection zone (e.g., 
Hindman et al.l 120041) . However, the agreement between the two analysis techniques cannot rule 
out that both techniques may be affected by the same or similar systematic effects. In this Let- 
ter, we report on a systematic center-to-limb variation in helioseismic travel times measured by 
the time-distance helioseismology technique, which was previously unnoticed but must be taken 
into account in the inference of the subsurface meridional flows. Other helioseismology tech- 
niques, such as the ring-diagram analysis, may also be affected by a similar systematic effect (Rick 
Bogart, private communication). We develop an empirical correction procedure by measuring the 
center-to-limb variation along the equatorial area during the periods when the solar rotation axis 
is perpendicular to the line-of-sight, i.e., when the solar B-angle is close to 0°. This correction 
scheme provides consistent results for the acoustic travel times in the north-south directions mea- 
sured from different HMI observables and AIA chromospheric intensity variations. We introduce 
our data analysis procedure and present results in §2, and discuss the results and their implications 
in 53. 
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2. Data Analysis and Results 
2.1. Data Analysis Tools 



Pesnell et al 



To fac ilitate analysis of the large amount of d ata from the Helioseismic and Magnetic Im- 
ager (HMI: IScherrer et al.ll2012l: ISchou et al.ll2012|) onboard Solar Dynamics Observatory (SDO; 



2012h . a time-distance helioseismology data-an alysis pipeline wa s developed and im 



plemented at the HMI-AIA Joint Science Operation Center (|Zhao et al.ll2012|) . Every 8 hours, the 
pipeline provides measurements of acoustic travel times, and generates maps of subsurface flow 
and wave-speed perturbations by inversion of the measured travel times, covering nearly the full- 
disk Sun with an area of 120° x 120°. 

For the pipeline processing, we select 25 overlapping areas on the solar disk and analyze 
the 8-hr sequences of solar oscillation data separately for each region. Acoustic travel-time maps 
are obtained for 1 1 selected wave travel distances, and are inverted to derive r naps of subsurfac e 
velocity and wave-speed perturbations from the surface to about 20 Mm in depth fehao et al.l2012|) . 
Then the results for individual regions are merged into nearly-full-disk maps covering 120° in both 
longitude and latitude, with a spatial sampling of 0?12 pixel^^ on a uniform longitude-latitude 
grid. The pipeline gi ves acoustic travel time measurements from two different fitting techniques 
dCouvidat et al.ll2012[) . and inversion results based on ray-path and Bom-approximation sensitivity 
kernels. Both measurement uncertainties for different distances and inversion error estimates for 



different inversion depths are given in IZhao et al.l (l2012h. In this Let ter, only the acoustic travel 
times obtained from Gabor-wavelet fitting (|Kosovichev & Duvalll997|) and inversion results based 
on the ray-path approximation kernels are presented. 

Although the pipeline is designed to analyze the HMI Dopplergrams, it can nevertheless be 
used to analyze other HMI observables that carry solar oscillation signals, e.g., continuum inten- 
sity, line-core intensity, and line-depth. T he full-disk data fro m the 1600 A and 1700 A channels 
of Atmospheric Imaging Assembly (AIA; Lemen et al.ll2012h onboard SDO can also be used for 
helioseismology studies with good accuracy ([Howe et al.ll201 lb. In this study, we compare results 
obtained from the four HMI observables and the AIA 1600 A data following the same analysis pro- 
cedure. This comparison helps us to identify the systematic center-to-limb variation and develop a 
correction method. 



2.2. Center-to-Limb Variation in Measured Travel Times 



The north-south (rns) and west-east (r„c) acoustic travel-time differences approximately rep- 
resent the north-south and west-east flow components, respectively, although a full inversion is 



required to determine more precisely these flows. We first show the measured travel time differ- 
ences and then present the inversion results. We choose a 10-day period of December 1 through 
December 10, 2010, when the solar B-angle between the equator and the ecliptic is close to 0° to 
avoid complications caused by leakage of the solar rotation signal into the meridional flow mea- 
surements. 

The upper panels of Figure [T] show the nearly-full-disk map of t^s, averaged over the 10-day 
period, measured from the HMI Dopplergram and continuum intensity data for an acoustic travel 
distance of 1?08 — 1?38. The general pattern of positive r^g in the southem hemisphere and negative 
Tns in the northern hemisphere is usually thought to be caused by the interior poleward meridional 
flows. However, the apparent differences between the magnitude of Tns obtained from the Doppler 
data and that obtained from the continuum intensity data indicate that there are additional system- 
atic variations. The lower panels show the averaged r^e maps measured from the same observables 
after the latitude-dependent travel times caused by the differential rotation, obtained by averaging 
the measurements of all longitudes, are subtracted. One would expect r„c be relatively flat along 
same latitudes because the solar rotation does not vary significantly with longitude, but both r^e 
maps show systematic travel-time variations along the same latitudes, positive in the eastern hemi- 
sphere and negative in the western hemisphere. The longitudinal variation is quite significant for 
the measurements from the intensity data (Figure [TJl) and small but not negligible for the measure- 
ments from the Dopplergram data (Figure [1]:). 

To more quantitatively illustrate this systematic variation, we average Tns in a 20°-wide band 
along the central meridian as a function of latitude, and display the averaged curves in the top 
panels of Figured for three selected measurement distances and for the four HMI observables: 
Doppler velocity, continuum intensity, line-core intensity, and line-depth. In the middle-row pan- 
els, we show the corresponding t„c curves obtained by averaging over a 20°-wide band along the 
equator as a function of longitude (hereafter, longitude is relative to the central meridian). If there 
were no systematic center-to-limb variation in the measured acoustic travel times, t^s obtained 
from the different observables would agree with each other, and r^e would remain flat for all ob- 
servables. But clearly, the measurements are not as expected. Among all these measurements, the 
Tns curves show not only different magnitude of travel-time shifts but sometimes also opposite vari- 
ation trends. For r„c, it is found that the measurements from the Dopplergrams have the smallest 
systematic variations and are similar to the line-core intensity measurements at larger distances. 
However, the travel times measured from the continuum intensity and line-depth data show not 
just substantial travel-time shifts for all measurement distances but also opposite center-to-limb 
variations. It is quite clear that the travel-time variations along the equator are not caused by solar 
flow but represent a systematic center-to-limb variation. If we treat the longitude-dependent r^g 
variations along the equatorial area as the systematic center-to-limb variations, and subtract these 
variations from the latitude-dependent Tns measured using the corresponding observables, we get 



the residual curves shown in the bottom-row panels of Figure [21 It is remarkable that for all the 
measurement distances, the results from all four HMI observables are in reasonable agreement. 
This suggests that the residual travel times correspond to the subsurface meridional flow signals. 

It was de monstrated that th e SDO/AIA 1600 A data are also suitable to perform helioseismol- 



ogy analysis (IHowe et al.ll201 1|) . Figure [3] shows a comparison of the results obtained from the 
HMI Dopplergrams and the AIA 1600 A data for one selected measurement distance, 3? 12 — 3?84. 
Similarly, Tns measured from the AIA 1600 A data also differ from that measured from the HMI 
Dopplergram. The r„c measured from the AIA data also show systematic variations, though its 
longitudinal trend is opposite to that measured from HMI continuum intensity and the variation 
magnitude is much smaller. The residual, obtained after subtracting the longitude-dependent Twe 
from the latitude-dependent Tns for AIA 1600 A, is quite similar to that obtained from the HMI 
Dopplergrams and other HMI observables. This result is particularly remarkable because HMI and 
AIA are different instruments observing in different spectral lines: AIA in 1600 A and HMI in 
Fe I 6173 A, which are formed at different heights in the solar atmosphere. 



2.3. Effect on Meridional Flow Inversions 

The systematic center-to-limb travel-time variations would have an apparent effect on the 
inference of subsurface meridional flows obtained by inversion. Here, we investigate how the 
inferred meridional flow speed changes after removal of this systematic effect. 

We employ the t ime-distance hel ioseismology pipeline code developed for inversion of the 



acoustic travel times (IZhao et al.ll2012h . The inversion results are shown in Figure |4] The middle 
row of Figure |4]presents inversion results obtained from travel-time measurements shown in Fig- 
ure |2l but with an extension of 15° closer to both poles and both limbs to better illustrate that the 
removal of the systematic effect also helps to improve inversions in high-latitude areas. The top 
and bottom rows of Figured show results obtained from June 1 - 10, 2010 and June 1 - 10, 201 1, 
when the solar B-angle were also close to 0°. It can be seen from left columns of Figure |4] that 
above ~ 55° latitude, the inferred meridional flow drops in speed and becomes equator-ward for 
the depth of - 1 Mm, but remains pole-ward deeper than 3 Mm. This is a suspicious behavior 
for the meridional flows. We then antisymmetrize the east-west direction velocity caused by the 
center-to-limb variations (middle column of Figure |4j, and subtract it from the inferred meridional 
velocity. The residual subsurface meridional flows show much more consistent behaviors at dif- 
ferent depths in high-latitude areas (right column of Figure IH). From Figure HI one can also find 
that the center-to-limb variation measured from equatorial area also slightly change with time, and 
why there is such a change is not understood and worth further studies. 
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The removal of the center-to-Hmb variation in the measured travel times also decreases the 
inferred flow speed by nearly 10 m s~^. This indicates that the subsurfac e meridional flows derived 
from the previous time-distance studies (e.g. Jzhao & Kosovichevll2004|) might have overestimated 



the flow speed by a large fraction. Figure |5] shows that after removal of the systematic effect, 
the meri dional flow speed is closer in m agnitude to the results obtained by the magnetic feature 



tracking (|Hathaway & Rightmirell2010t) and surface Doppler measurements (lUlrichll2010l) . Note 
that the results from the magnetic feature tracking and from the MWO Doppler observations shown 
in Figure|5]are both averaged over a 6-month period with our analyzed 10-day period in the middle. 
The difference of analysis period may result in some differences seen in the figure. The meridional 
flows obtained from this study are displayed in Figure |5] after a 2° spatial averaging to remove the 
strong fluctuations caused by supergranular flows. 



3. Discussion 

The analysis of acoustic travel times obtained from different observables of the HMI and AIA 
instruments on SDO by the time-distance helioseismology technique has revealed a systematic 
center-to-limb variation. The systematic variation is different for different observables, and range 
from ~2 sec for the HMI Dopplergram to ~10 sec for the continuum intensity measurements. For 
an accurate determination of subsurface meridional flows, and also for more accurate inference of 
full-disk subsurface flow fields, this systematic effect should be removed. We have developed an 
empirical correction procedure by removing the systematic variation measured along the equatorial 
area during the period when solar B-angle is close to 0°. This correction reconciles the latitude- 
dependent Tns measured from different observables and reduces the inferred meridional flow speed 
by about 10 m s^. 

It is not quite clear what causes this systematic center-to-limb effect. However, this effect is 
unlikely caused by instrumental or data calibration, because it is observed in the data from two 
different instruments, HMI and AIA, and it exists in different observables of HMI. We also rule 
out the following factors as causes of this effect: finite speed of light, contribution of horizontal 
wave component, and foreshortening effect. Due to the finite speed of light, acoustic wave signals 
observed at high latitude and observed near the equator are not simultaneous, but this effect is neg- 
ligible given the srn all measurement distances used in this study. The horizontal wave component 
(INigam e t al. 2007; is used to explain the center-to-limb variation in mean travel times observed 



by iDuvalL (^2003i) . but this effect is not significant in travel-time differences for waves traveling in 



opposite directions and in small distances as used in this study. To check foreshortening effects, 
we mimicked the high-latitude data by reducing the spatial resolution of the data observed at lower 
latitudes, and our measurements did not show systematic changes in the measured travel times 
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similar to the center-to-limb variations reported here. 

It is more likely that this effect has a physical origin related to properties of solar acoustic 
waves in the solar atmosphere or response of the spectral lines to the solar oscillations. It is well 
known that for a given spectral line the Sun is observed at the same optical depth but not at the 
same geometrical height. The observed height gradually increases with distance from the disk 
center. Thus it is possible that the acoustic waves traveling in opposite directions and different 
height will give different measured travel times due to the subsurface location of the wave source 



and the waves evanescent behavior above the photosphere (|Nagashima et al.ll2009i) . However, this 
cannot explain why different observables give different magnitude and trend of the center-to-limb 
variations. Another interesting fact is that the largest systematic effect exists in the measurements 
using HMI continuum that forms at the lowest height in the atmosphere. As one moves up to 
the height where Doppler velocity is measured, the systematic effect is reduced. As the measure- 
ments move further up to where AIA 1600 A is formed, the effect reverses sign. This may give 
us some indication of that this center-to-limb variation is related to the line formation height. It is 
also possible that this systematic effect is related to differences in the acoustic power distributions , 
line-a symmetry of solar modes f Puvall et al.lll993|) . and the correlated noise effect (iNigam et al 



1998i) . The cause of this systematic effect is worth further studies, and may be resolved by numer- 



ical modeling of solar oscillations including the spectral line-formation simulation and spherical 
geometry of the Sun. It is worth attention that when inferring meridional flow velocity from sur- 
face Doppler measurements, a sy stematic center-to-limb velocity profile also needs to be removed 



(lUlrich et al.lll988l: IUlrichll2010|) . A similar limb-shift effect very likely exists in the HMI Dopp- 
lergrams that are used in our time-distance measurements, and it is not clear whether this effect 
accounts for some of the systematic variations in our measured acoustic travel times. This is worth 
further studies. 

Although the cause of this systematic center-to-limb variation is not well understood, it is 
demonstrated that the empirical correction procedure can improve the inferred subsurface merid- 
ional flows. Figures |2] - ID demonstrate that subtracting the longitude-dependent r„o from the 
latitude-dependent Tns, measured from same observables, is an effective way to remove this sys- 
tematic effect. This correction procedure helps to reconcile the Tns measured from different ob- 
servables, and also helps to remove the inconsistent behaviors of meridional flows in high-latitude 
areas. As an effect, the newly obtained subsurface meridional flows at shallow depths are approxi- 
mately 10 m s^^ slower than the speed previously derived following a similar analysis procedure. 
This systematic center-to-limb variation in measured acoustic travel times has an important im- 
plication in the long-searched deep equator-ward meridional flows, and this is currently under 
investigation. 
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Fig. 1. — Travel-time maps (in unit of sec), covering 120° in both longitude and latitude, for r^s 
obtained from (a) HMI Dopplergrams and from (b) HMI continuum intensity, and for t„c obtained 
from (c) Dopplergrams and from (d) continuum intensity after a longitudinally averaged profile, 
representing the differential rotation, is removed through all longitudes. Note that color scales are 
not the same in all panels. 
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Fig. 2. — Top: Averaged curves of latitude-dependent Tns, obtained from different HMI observables 
and for different measurement distances. Middle: Averaged curves of longitude-dependent r„e, 
obtained from different observables and for different measurement distances. Bottom: Differences 
of Tns and r„e. Note that the vertical scales for the upper two rows are different from those for the 
bottom row. 
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Fig. 3. — (a) Averaged latitude-dependent r^s, obtained from HMI Dopplergrams and from AIA 
1600 A data for a measurement distance of 3? 12 — 3? 84. (b) Averaged longitude-dependent r„e for 
these two observables. (c) Differences of t^s and r„e. 
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Fig. 4. — Results obtained from the periods of June 1 - 10, 2010 (top row), December 1-10, 
2010 (middle row), and June 1 - 10, 201 1 (bottom row), (a) Meridional flow velocity for different 
depths, obtained by inversions of acoustic travel times measured from the HMI Dopplergram and 
without removal of the systematic center-to-limb variations, (b) Antisymmetrized east-west veloc- 
ity obtained by inversion, representing the center-to-limb variations, (c) Meridional flow velocities 
after removal of the systematic center-to-limb variations, (d) - (f) and (g) - (i) are the same as (a) 
- (c) but for different time periods. 
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Fig. 5. — Comparison of th e meridional flow speed obtained by magnetic fe ature tracking 
(IHathaway & Rightmirell2010h . from surface Doppler measurements (|Ulrichll2010|) . and obtained 
in this study for a depth of - 1 Mm before and after the removal of the systematic effect. Our 
results are from a 10-day period, December 1 - 10, 2010, and are displayed after spatial averaging. 
The results from magnetic feature tracking and from Doppler measurement are averaged over a 
6-month period with our analyzed period in the middle. 



